Measurement of the Slope Parameter for the rj — > 3ty° Decay in the pp — > ppr] Reaction 
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The CELSIUS /WAS A setup is used to measure the 3n° decay of n mesons produced in pp inter- 
actions with beam kinetic energies of 1.36 and 1.45 GeV. The efficiency-corrected Dalitz plot and 
density distributions for this decay are shown, together with a fit of the quadratic slope parameter 
a yielding a = —0.026 ± O.OlO(stat) ± O.OlO(syst). This value is compared to recent experimental 
results and theoretical predictions. 

PACS numbers: 13.25.-k, 14.40.Aq, 12.39.Fe, 13.75.-n 
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I. INTRODUCTION 



Isospin violating decays into three pions are two of the 
most probable n meson decay modes: BR(rj — > 37r°) = 
32.5% and BR{r/ -» ttV+tt") = 22.6% [1]. It appears 
that contribution of electromagnetic processes is sup- 
pressed 0, H[ and the decays are driven by an isospin 
violating term of the QCD Lagrangian proportional to 
TOrf — m u . The partial width of the r\ — > 7r°7r + 7r~ de- 
cay calculated using current algebra (PCAC) is 66 eV 
[|J, much below the experimental value 294±16 eV [jj. 
Second order contributions in the low energy expansion 
of the QCD Lagrangian (Chiral Perturbation Theory - 
CHPT) were calculated by Gasser and Leutwyler increas- 
ing the result to 160 eV [5J. The big change with respect 
to the first order calculations implies the importance of 
7T7T interaction in the final state. An elegant method of 
including the interaction up to higher orders is provided 
by dispersion relations. There are two calculations using 
this technique @, 0] but employing different formalism. 
They lead consistently to an enhancement of the decay 
rate by about 14%. A free parameter in this approach 
is the value of the so called subtraction point and it was 
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constrained using CHPT calculations. 

The decay width can be expressed in the factorized 
form: 



F 



Qd 

Q 



f (1) 



where the dependence on the rrid 
in the Q term: 
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Normally the calculations of T are performed using the 
Dashen theorem Q where Q — Qd = 24.1. Since T is 
sensitive to the exact value of Q, it was suggested that 
the decay might provide precise constrain for the light 
quark mass ratios @. Namely Q determines the major 
axis of the ellipse in the m u /m Sl md/m s plane. One 
important prerequisite is the reliability of the T calcula- 
tions. The calculations can be tested by comparing the 
ratio T(r) — > 37r°)/r(?7 — > 7r + 7r~7r°) and the density dis- 
tributions in the respective Dalitz plots with experiment. 
Recently Q was derived from preliminary KLOE data on 
r] — > 7r + 7T~7r° decay [lfjl [j"lj|. by determining the subtrac- 
tion point within the dispersion relation approach of 0], 
yielding a Q value 22.1 [121 ]. 

This letter describes an exclusive measurement of kine- 
matical distributions of the ir° 's from the r\ — > 3tt° chan- 
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TABLE I: Experimental and theoretical results for the slope 
parameter a 



a 


Comment 


Ref. 


-0.052 ± 0.017(stat) ± O.OlO(syst) 


Exp (CBarrel) 


M 


-0.031 ± 0.004 


Exp (CBall) 


[19] 


-0.013 ± 0.004(stat) ± 0.005(syst) 


Exp (KLOE) 


[10] 





PCAC 


1 


+0.015 


CHPT,Hoop 


[5, 20] 


-0.007.. .-0.014 


CHPT+dispersive 


w 


-0.007 


UCHPT 


[13] 


-0.031 ± 0.003 


UCHPT/fit 


[21] 



nel. The phase space of the decay can be represented by 
a Dalitz plot employing symmetrized variables as shown 
in Fig. QJleft). The amplitude can then be expressed in 
terms of polar coordinates p and 0, where due to sym- 
metry (f> is restricted to < cj) < 60°. Instead of p a 
normalized dimcnsionless variable < z < 1 is used: 
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where Ti is kinetic energy of a 7r° in the r\ rest frame and 
< T > is the average kinetic energy. 

The slope parameter, a, is defined as a coefficient in 
the leading term of the decay amplitude expansion in z: 



A(z,4>) \ 2 = c a (l + 2az). 



(4) 



The situation for the a measurement is unsettled with 
inconsistent experimental results and disagreement with 
CHPT predictions (Table 0}. In lowest order CHPT 
a = 0. The deviation is caused by the 7r°7r° interactions. 
All calculations beyond one loop predict a negative value 
of a. Dispersion relations H an d U(3) chiral effective 
field theory with rescattering jl3| give consistently a low 
value for the |a| of about 0.007. Most recent calcula- 
tions with U(3) effective Lagrangian and Bethe-Salpeter 
equations include more general rescattering graphs and 
parameters are fitted to the scattering data for the pseu- 
doscalar mesons [l4| • The calculations are able to explain 
the large experimental values for \a\. However in that ap- 
proach a change of the Q value can be compensated by 
other fit parameters and could even be equal to Qd- 

A special feature of the 77 — > 37r° decay is that the 
physical region extends below threshold of the 7r°7r° — > 
7r + 7r~ reaction. Therefore a cusp structure should be 
visible in invariant mass distribution of two ir° around 
2m^± (dashed lines in Fig. DJleft)), in analogy to the 
recent observation in K + — ► 7r + 7r°7r° decay [HdCLlI]. 
The exact parameters of the cusp in the K + — > 7r + 7r°7r° 
decay provide one of most precise determinations of the 
7T7T scattering length (ci2 — ao)- 

The experimental determination of the a requires that 
background contribution to the i] production process and 
other systematic errors are well under control. All recent 
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FIG. 1: (Left) Symmetrized Dalitz plot for r\ — » 3n° decay. 
(Right) radial density distribution (right) for the MC (shaded 
area) and experimental data (closed symbols). 



experiments on the decay rj — » 37r° rely solely on pho- 
ton detection with close to 47r electromagnetic calorime- 
ters. The Crystal Barrel experiment used the final state 
with IO7 (pp — > t]2tt°) [18j, Crystal Ball used the reac- 
tion n~p — > nrj where only the 67 from the decay were 
measured [ljjj]. KLOE [l(| used the decay 4> — > 777 as a 
source of the 77 mesons leading to 77 in the final state. 
The a values differ by more than three standard devia- 
tions between the experiments. We have performed an 
independent measurement using the reaction pp — > pprj 
[22I ]. The main advantage is that the excellent missing 
mass resolution of the forward going pp pair can be used 
to tag the process. 



II. MEASUREMENT AND DATA REDUCTION 

The analysis is based on data taken with the WASA 
detector at CELSIUS [23| using a pellet target system 
that provides small (30 /im diameter) hydrogen pellets 
that cross the proton beam of nominal 1.36 GeV and 
1.45 GeV kinetic energy (the corresponding center of 
mass excess energies are 41 MeV and 75 MeV respec- 
tively). The WASA detector system comprises a mul- 
tilayered forward detector (FD) for the measurement of 
charged particles emerging in the scattering angle range 
of 2.5°-18°, and a central deteector (CD) composed of 
an electromagnetic calorimeter of 1012 CsI(Na) crystals 
and a drift chamber/solenoid combination for the mea- 
surement of particles from rj decays in the angular range 
of 20°-140°. 

The basic criteria for the selection of the final 77 — > 3ir° 
state are (i) two protons detected in the FD and (ii) six 
7 hit clusters (from 7r° — > 77) detected in the CD as 
neutrals with a minimum of 20 MeV deposited energy. 
The experimental proton-proton missing mass distribu- 
tion can be compared to Monte Carlo (MC) simulated 
data for the reaction channels pp — > ppn c 'ir c 'tt° , pp — > ppr], 
an d pp — > p7J7r 7r°. The procedure is described in detail 
in [22|, the result for 1.36 GeV kinetic energy is shown 
in Fig. O The good agreement in shape that is obtained 
also for other kinematical proton variables and of the 
absolute cross sections for the 77 channel with existing 
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FIG. 2: Experimental pp missing mass distribution for recon- 
structed ppir°ir°n° final states with a fit of the MC distribu- 
tions for 27r°, 3ir° and r\ production (from [22|]). The vertical 
lines indicate the r\ selection. 



data, demonstrates that the detector and its efficiency 
are well under control. The missing mass resolution of 
the FD (5 MeV/c 2 FWHM at 1.36 GeV) therefore allows 
selection of the pp — > pprj reaction (the condition 535 
MeV/c 2 < MM pp <560 MeV/c 2 was applied) such that 
the remaining background, mostly from direct pp — > pp3n 
reaction, is about 5%. 

All possible combinations (15) of the 6 reconstructed 
gammas to form 3 7r° — > 77 pairs are sorted by means 
of a parameter x 2 = f=i(-^-^i — rn^o)" 2 ''/erf, where IMj 
is the invariant mass of the i-th pair of a combination, 
and Ci the resolution. At most two combinations with 
the lowest \ 2 are selected for a kinematical fit of the full 
event with 8 constraints: four-momentum conservation, 
the 7r° — * 77 constraints and the rj — > 37r° constraint. 
A cut on the \ 2 of the most probable combination is 
applied to further suppress background and to increase 
the combinatorial purity. Three different data sets were 
individually analyzed, the combined result is based on 
75000 events in the Dalitz plot after all cuts. 



III. DALITZ PLOT AND SLOPE a 

The reconstructed experimental z distribution is shown 
in Fig. QJright) together with MC simulation including 
pp — > pprj and pp — > pp7r 7r 7r° reactions. The detector 
response function for the reconstructed z was determined 
from MC: it has a gaussian distribution with a standard 
deviation (RMS) 0.067. The resolution in the 7r°7r° in- 
variant mass (-Mtttt) is approximately constant and equal 
to 6 MeV/c 2 (RMS). 

The acceptance corrected |A| 2 dependence on the z 
variable is obtained by dividing the measured distribu- 
tion by the MC prediction with a set to 0. A fit of 
the function Q is performed to extract the parameter 
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FIG. 3: Extracted experimental |A| 2 /cn dependence on the z 
variable for a subset of the T = 1.36 GeV (left) and for all 1.45 
GeV (right) data together with a fit of the slope parameter 
a. Only statistical errors are shown. 



TABLE II: Statistical errors of a fitted for data subsets ob- 
tained by parameter variation for selection criteria and recon- 
struction procedure. 



Condition 


a± 


O stat 


events 


All data, 0.1< 2 < 0.9 


-0.026 


± 


0.010 


74700 


All data, 0.0< 2 < 0.9 


-0.014 


± 


0.009 


74700 


Analysis of subsets: 










Subset I T = 1.36 GeV 


-0.023 


± 


0.018 


20700 


Subset II T = 1.36 GeV 


-0.041 


± 


0.018 


18000 


Subset T = 1.45 GeV 


-0.019 


± 


0.015 


36000 


Subset I, variation of MM VV cut: 










0.530 < MM PP < 0.575 


-0.019 


± 


0.018 


21300 


0.535 < MM PP < 0.560 


-0.023 


± 


0.018 


20700 


0.540 < MMpp < 0.555 


-0.025 


± 


0.018 


18800 


Subset I, variation of Xkinfit cn ^ : 










Xlinfit < 999 


-0.023 


± 


0.015 


27900 


Xkinfit < 50 


-0.021 


± 


0.017 


23600 


Xkinfit < 30 


-0.023 


± 


0.018 


20700 


Xkinfit ^ 15 


-0.021 


± 


0.020 


15200 


Xkinfit — 15; Xsecond/X. first — l - ^ 


-0.038 


± 


0.024 


10800 


Subset I, alternative parameterization 










of kin. fit 


-0.023 


± 


0.019 


17000 
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slope a= -0.026 ± 0.010 (stat) 
based on 75000 events 
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FIG. 4: Extracted experimental |yl| 2 /co dependence on the z 
variable for all data together with a linear fit to of the slope 
parameter a. Only statistical errors are shown. 
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FIG. 5: Ratio of the M,^ distribution to the MC predictions 
with a = 0: points - data; line - CHPT calculations from 
reference Hal. 



a. The MC data are normalized to obtain cq — 1. For 
the final result a = -0.026 ± 0.01 (stat) ± O.Ol(sysi) the 
first bin was excluded due to large systematic uncertainty 
and the last bin due to low statistics. The estimate of the 
systematic errors is obtained by variation of all essential 
cuts applied in the reconstruction: \ 2 f° r the kinematic 
fit, the combinatorial purity of the sample, the missing 
mass and the z range. For each subsample a fit to a was 
performed: an example is shown in Fig. [3] and the sum- 
mary is given in Table [Til The overall systematical error 
was obtained by comparing central values of the fits for 
the subsamples and it takes into account the influence of 



the first bin for the a value [241 ]. 

We have also performed a search for a cusp structure 
in the invariant mass distribution of the two 7r° (Fig. [5]). 
An eventual observation would require at least two orders 
of magnitude larger data sample. 

The presented result for a is limited by the available 
statistics which determines also attainable systematical 
accuracy (due to the size of the subsamples used for the 
tests). Within the errors it is compatible both with the 
result from Crystal Ball [lj|, based on 10 6 events, and 
with the preliminary KLOE result [H| listed in Table [TJ 
The reconstruction methods developed for the purpose of 
this experiment are now being used in analysis of recently 
collected data from the followup experiment with WASA 
detector located at COSY ring (FZ Juelich) (25| with 
much larger statistics. 
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